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УДК 539.3 
Колебания неоднородного пороупругого слоя с пустыми порами! 


А. О. Ватульян, Д. С. Шведов 
(Донской государственный технический университет) 


В рамках плоской деформации рассмотрена задача об установившихся колебаниях неоднородного по тол- 
щине изотропного пористого упругого слоя с пустыми порами. Нижняя грань слоя сцеплена с абсолютно 
жёстким основанием, колебания вызываются поверхностной нагрузкой на верхней грани. При известных 
законах неоднородности для аналогов параметров Ляме — положительных функциях вертикальной коорди- 
наты — колебания слоя описываются системой трёх уравнений в частных производных второго порядка с 
переменными коэффициентами. С помощью преобразования Фурье и обращения некоторых операторов за- 
дача сведена к системе трёх интегральных уравнений Фредгольма второго рода с непрерывными ядрами. 
Предложен численный метод нахождения трансформант на основе метода коллокаций. С помощью числен- 
ного обращения преобразования Фурье осуществлён расчёт вектора смещения и функции относительного 
объёма. 

Ключевые слова: пороупругость, система интегральных уравнений Фредгольма второго рода, функция 
относительного объёма. 


Введение. Теория пористых линейно упругих материалов с пустотами первоначально была раз- 
работана Коуином и Нунциато [1]. Она предназначена для моделирования упругих материалов, 
содержащих распределение малых пор (рис. 1). Главной особенностью данной теории является 
введение в определяющие соотношения новой переменной, характеризующей относительный 
объём пор, который берётся в качестве независимой кинематической переменной. Включение но- 
вой переменной требует дополнительных сил для обеспечения равновесия объёма пор. Если объ- 
ём пор обращается в ноль, поведение материала описывается классической теорией упругости. 


Порода материала 


Рис. 1. Упругая среда с распределёнными порами 


Теория Коуина—Нунциато хорошо описывает поведение композитных, керамических, пе- 
нистых и пористых материалов. При этом ряд функционально-градиентных материалов отличает- 
ся тем, что распределение физических характеристик неравномерно, что приводит к их зависимо- 
сти от координат. 

Настоящая работа посвящена исследованию колебаний неоднородного по толщине поро- 
упругого слоя с пустыми порами. 





1 Работа выполнена в рамках инициативной НИР. 
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Постановка задачи. Рассмотрим двумерную задачу об установившихся колебаниях с частотой © 
неоднородного изотропного пористого упругого слоя 0 < у< А, |х| < ©. Нижняя грань слоя гори- 
зонтальна и сцеплена с абсолютно жёстким основанием. 

—ГОЕ 


Необходимо найти поля перемещений в слое и(х, у) ={и,и}е и функцию относитель- 


ного объёма Ф(х,у,Ё) =Ф(х,у)е"" при известных законах неоднородности под действием 
нагрузок, приложенных к верхней границе слоя. 
Считая режим колебаний установившимся, отделим временной множитель е`^*. Будем 


рассматривать амплитудные значения функций. 
Уравнения движения имеют вид [2]: 


2 
би х +07, +Р®"И =0, 




















1 
бк +0,, +РО№ =0. (1) 
Компоненты тензора напряжений находятся из следующих определяющих соотношений [2]: 
9х ай 2 ' 
мт =, + (1-2с у, + НФ, 
б (2) 
У = 2 , , 
т = (1-2с? и, +и, + НФ, 
бу =Н(м, и} 
где Н = В ‚1-26? = Л : 
А+2Н А+2Н 
Полевые характеристики задачи должны быть подчинены следующему уравнению [2]: 
2 2 : 2 
Я Сы и В ди ‚ду 0 (3) 
0х ду а ах ду 
Соответствующие граничные условия для слоя могут быть записаны в виде: 
и(х,0) =0, 9, (х,й)=т(х), 
у(х,0) =0, о„(х,В)=ов(х), (4) 
5 (х,0)=0, |2 (ж,в)=0, 
бу бу 


где Л =^(у), и= н(у) — положительные функции координаты у, являющиеся аналогами пара- 
метров Ляме, р — массовая плотность материала; а, В, Е, ®, К — положительные физиче- 
ские параметры, связанные с пористостью [1]. 
Система интегральных уравнений Фредгольма 2 рода. Отметим, что в представленной по- 
становке колебания описываются системой уравнений в частных производных второго порядка с 
переменными коэффициентами. Наиболее эффективным средством анализа колебаний в слое яв- 
ляется использование преобразования Фурье по переменной х и сведение к системам обыкновен- 
ных дифференциальных уравнений с переменными коэффициентами [3]. 
Функция 
(= [7 (%)е®ах 
называется прямым преобразованием Фурье Ё (х). 
В свою очередь, функция 
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называется обратным преобразованием Фурье и с физической точки зрения это выражение 
представляет собой интегральную запись принципа суперпозиции гармоник с комплексной 
амплитудой и волновым числом $ с непрерывным спектром. 

Системы (1), (2) и уравнение (3) в образах Фурье выглядят следующим образом: 


Ри а - ы 
—56,. +— 6, +рО?й =0, 
хх бу ху р 











(5) 
5 9:5 ОУ =0 
О у +В и =0, 
6, =-ва (А+ 26) + (и 2с*(А +2) уч НФ (А+), 
ы ыы Не 
0 = (1-2с*)(^+25)(-5)й и + №Ф (+24), (6) 
б, =Н ба [5й ЕЙ 
г ау ау | 
* 7 2, _ 
об В В Чт а |0. (7) 
ау а а\ ау 
2 
Если из (5)—(7) выразить функции г б ИИ Ф, исключить б„, а также 





б у 

ау ”'ау ”ду ‘ду 'д’ 

а: - ы 

ввести о = Ф , получается каноническая система обыкновенных дифференциальных уравнений 
У 


1-го порядка с переменными коэффициентами: 














Чел в Ре и(А+ы) о 
6. =5 б ©? — 45? й +525, 
ау ” л+2н ” р А-+2Н не 
я 5, =0,„, -РОМ, 
Ч 1-16 +5 
Е (8) 
а - й- $ 
а  М№+2и” л+ | 
Ч ф-$, 
ау 
— [ 2%. — 
ее Аи оо ВЕ Ил РАНЫ ВЫ 
ау а(^+2н) ” ал+2 а а 


Граничные условия в образах Фурье имеют вид: 
0 (5,0) =0, б„ (5,№)=1(5), 





ху 

о _ ‚й)=6(5), (9) 
,0)=0, ‚В )= 
(5.0) -0, [9(51)-0 
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Краевая задача для системы функций (8), (9) после введения безразмерных переменных и пара- 


метров приведена к следующей безразмерной форме: 





2-, “14 (уе [0,в]52=[01], 











в’ ау Ва 
9, (избы, ©, (избы), И (5,2) = 
(в, 2) 782), (5,2) =$(5,2), 
(2-е), 
ке рада 
т(5) 








Й 


После замены для получения вещественных коэффициентов в вырожденной системе (& =0) 


[И =(, Ю =0,, 9 =0,, 


она выглядит следующим образом: 














а х н, (^, +1) 921? 

=-№5 ——\ д И —- 2654. Н.Ф, 
42 © Л, +25, © + Л +24, И Вр 
а р@?8? 
— О, = 150, - у, 
42 ©, ©, № 
а 1 
ЕО ВЫ, 
|= НбНТ м 
а 1 К 
0 Но 
2 М -+2ы © + А, * ‚ФР 
Ч ф-у 
[е7. ' 
а р? 

В [] 





@ а(л, +24) а 


Граничные условия: 


а Л. +24 


И, (5,0) =0, |0, (5,1) = Ло (5), 
У (5,0) =0, 10, (5,1)=0. ($), 
4(5,0)=0. |4(5,1)=0. 


В [0] |” №0 +рКо? -Е 25? 
1 


(10) 


(11) 


р? 
РН 


(12) 


В общем случае построение решения канонической системы дифференциальных уравнений (11) с 
переменными коэффициентами в аналитической форме невозможно, поэтому осуществляется её 
сведение к системе интегральных уравнений Фредгольма второго рода. Интегрируя систему (11), 
учитывая граничные условия (12) и осуществляя элементарные преобразования, получим систему 
уравнений Фредгольма 2 рода с непрерывными ядрами (отметим, что все неизвестные функции 


зависят от параметра преобразования Фурье): 
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©, = [(6,0, (п) + СЦ, (п) Е С.Ф (п))ап = Пт, (1 — 2) +0, 


и, =-| (6.0, (п) + 6, (п) + 6:$(п))ап+ тии (2), (13) 





Ф= Г б,О, (п) + 6.0, (п) + 6. (п))ап, 


И 
Ядра системы имеют вид: 
Л, 021? 1 
1 Е 2! 
А, +2н, М Л +24, 
№ (А, +8.) р? |, 95° Л, 
А, +2, №] о А+2н 








2' 


22 
6, (2,п) =-252° В.Н, К, - м НК,, 


0 











1 
6. (2,п) =-1$ —— К, +15-—— К), 
« (2,1) А, +2н, з" ар“ 
№ ( +щ) р9?/? 22 
С. (2,1) =| 45202 о К. +528 1 ы 
в (Е,п) = | Л, +2, № - „+24, “ 


6. (2,п) = 2856. Н.К. - И5Н,К, 
2 Ее 
ЕН 

















а(^, +2ы, ) / 
ВА? Н 1-К; - К 
6,(2, 25 и © 
в (2,п) а Л, +24, ее 1 
р 2 _ 2 
б. (2,П) = и : + и, ке 


где 
Я — ы; = (п), 
К, (2,п) = тах {2,1} -2 
1 
К, (2,1)= [| рас, 


пах{2 1} 
пи1{ 71} 1 


Кз (2,п) = | — ас, 
0 ы, 
К. . о п}, 
9 (у) аъ, 


и} 


м 


о 
|“ 


9 (у) аъ, 


пах{2,п} 


1 2 2 
о №2О + — ив м к, 
П 


ре |1. ас, 


0 Н, 


1 р 2 
.- |" №0 + ро? —& ее ВИН ат 
0 











а 
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Численный метод решения системы интегральных уравнений и обращения преобра- 
зования Фурье. Осуществим дискретизацию системы интегральных уравнений Фредгольма 2 
рода (13), используя квадратурную формулу трапеций [4, 5] и метод коллокаций. Придём к сле- 
дующей алгебраической системе: 





и. = 5-Е Цеци, 
п-1 
А,=А, = =, Д-АЕКЕА, АХ, 
Я = УА, (610. + С, ке СФ, ) = 0, > 1И5Т, (1 И 74 ), (14) 
1 =1 
И, + >^, (6.,0,, +6, +быФ, ) = и, 
= 
ф,- >А (6,0, + 6:0; +б,Ф,)=0. 
Л= 





Значения остальных функций в соответствующих узловых точках О;‚,И,щ, находятся из инте- 
гральных представлений: 


РО л (А, +Н 021? 
О, (2) = То Г 5 ИХ 2 ©, + и "В р ол + 2н :) Е 
2 1 








|4 285р.Н‚ф | ас, 
1 №о 


т А 
+ 5—0 нев 
> +2, ©, + +2 г 1 
г Вл? Вин, › ЮО + ро? -Е „>, ВА" 
2) = 25 ив 975+ Н : 
+=) РЗ а Л+2ы ' а а” 


0 





Таким образом, для любых значений параметра преобразования можно найти узловые значения 
переменных, и далее требуется построение обратного преобразования Фурье. Для численного 
осуществления обратного преобразования Фурье, т. е. для вычисления интеграла 


х) == | Е ($)е 4$, 


используем интерполяционную формулу третьей степени, являющуюся аналогом «правила трёх 
восьмых» Ньютона-Котеса [5, 6]. Оно получается при замене функции Ё на каждом из отрезков 
[о,з^, | ,[З№,,6А, |,... алгебраическим многочленом третьей степени, интерполирующим Е по че- 
тырём значениям. 

Возьмём четыре точки $,,5,.1,5,,›,5,.з И ВЫПОЛНИМ Интерполирование / по её значениям 
в этих точках: 


Е (5) = (5 Ее - 5:2) ($ — 53) (5 -5%)(5 - 54:2) ($ — 53) 











[в (228, (31, А: 
(55%) (5 — 54) (5 — 52) = ‚ (6-54) ($ - 51)($ — 5.2) + (5). 


28, -В, (-В, т ЗА, 21, - №, _ 
Умножение этого равенства на е`^ и интегрирование по отрезку Е + ЗА | с последующим 


суммированием по значениям К = 0, +3, -6,... приводит к представлению: 


е- | — ©.) = 0.) = , 
х) => [Е (5)е°% = (У. = ‚) >. Ее + о +1 О РА 23 + а > Е „273 + А (х) ) (15) 


К=—о К=Ь—2о К=Ь—2о 2п 
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где 


1 


-* (6.2), -З)-уЕ-т,)+ 


в(х)= . [авт |6 5.) Ее) 











+5, (6 (=) =) Е@=т) 5, (6, 1) (5, -2)(3 ый 


1,5 > 0, 
($) = 5,5 =0, 


0, $ <0. 


Коэффициенты а,,у,,б, имеют вид: 
9=х№,, 


ва, = 8? 6 + [© в 9? 05 36 + 03 536, 
В ‘у, =30“ -30” +3 [ 9. Ос 30 — 40-3 5тз6, 


в.15, -503 +3 (Е 9? т 36 + 40 с0$36, 


Стоит отметить, что численный эксперимент проводился для л = 30 узловых значений при точно- 
СТИ & = 0,001. 

Расчёт полей смещений и функции относительного объёма при некоторых законах не- 
однородности. В качестве примеров реализации предложенной схемы рассмотрим возрастаю- 


щие и убывающие функции, характеризующие законы неоднородности в случае сосредоточенной 
и распределённой нагрузки. 


Для сосредоточенной нагрузки находим: 
о, (х)=5(х); в, (5) = | б(х)е®ах =1, 


© 


т (х)=5(х); ть (5)= Гобетах и 


Для равномерно распределённой на отрезке нагрузки получим: 


в. (Х)=5.; в аа Вы 
© (х)=0, |х|>а 5 
= я: Е ыы 
% (х)=0, |х|>а 5 


где а — полудлина области приложения нагрузки. 


Будем рассматривать законы неоднородностей, для которых выполняются следующие ра- 
венства: 





1 1 
[^ (2) а =1, Гы (2) =1, р(2)=1. (16) 
0 0 
Согласно физическому смыслу параметры Ляме должны удовлетворять условию: 
Л, (2) 1 
0 < - <= \2е[0;1|. (17) 
2(^, (2) +1 (2)) 2 [ой 
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Рассмотрим некоторый закон, удовлетворяющий условиям (16), (17): 
А, (2)=2+0,5, в, (2) =1,5-2. 
Пусть А =0,5; а=0,01; В=10; «=0,01; &=15; К =0,01; © =100 , атакже 
(5) [5 =0, 
Ниже приведены графики полей перемещения и функции относительного объёма. Сплошной ли- 
нией обозначена действительная часть, пунктирной — мнимая часть. 



































/\ 
Ил 
А 
Ра 
Л 
То | р | 16 20 
х р. 
——— КВе(о) — — порах ——— Веб == неСР | 
0,1 = 
/\ 
0 
—0,1 
—0,2 
—0,3 
—0,4: 
—0,5 
—0,6 
—0,7 
| | | . | | | 
—20 _10 [и 10 20 





2 
[—— вефь —-— нофь | 





Рис. 2. Поля перемещений в слое и функция относительного объёма 


Заключение. Разработан метод решения краевой задачи для неоднородного пороупругого слоя, 
позволяющий находить поля внутри и на поверхности для широкого класса неоднородностей. 
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ОЗСТЕАТТОМ$ ОЕ 1МНОМОСЕМЕОЦ$ РОВОЕГА$ТТС ГАУЕВ \М/ТТН \ОТО$'! 


А. О. Уауапт, О. $. $Нуедоу 
(Роп Зае Тесптка! Упмег$Ку) 


Аз рагЕ оЁ {те р/апе хат, те ргоМет оп $еаау озс!айоп$ о! те твотодепеоиз гоидв-Ш/сКпез$ 15отор/С 
рогое/а$Нс 1ауег и уоГа$ 15 сопз/АЕгеа. ТПе [вуег Бойот 15 соппесеа ий те ретесйу поа Гоипаайоп, Те 
о5сайоп$ аге саизеа Бу {Те 5итТасе 5'ез$ оп {Те юр асе. ИпаЕг {те кпоит ппотодепейу (ви/5 Юг [ате ра- 
гатеег апа/од$ — роз#ме Кипсвоп$ оЕ те уегиса! соогтейе -— ве Гауег осайоп$ аге аезстреад Бу {пе учет оЁ 
{ргее рагИа/ АИТегепва! зесопа-огаег едиайоп$ ий уапа Бе соейс!етё. Изтд {пе Еоипег Еап5Югт апа 5оте 
ассез5 ${ейетепё, {Те ргоМет 1/5 гедисе4 {о {пе зубет о!Ё {пгее Егеабо/т и{едга/ едиавоп5 ог {Те 5есопа кта 
ИИЁВ сопйпиоиз Кегпе. Тйе питепса! {есбтдие 15 ргорозеа Гог аепйутпа е ЕапзЮгт$ Бу те со!осаНоп 
тепоа. Тпе @зр/асетепЕ иесюог апа {Те ге/айуе иоите аге сасшаеа изтд питепйса! пуегзе Роипег вапЮгт. 
Кеуигогаб: рогоеа5Исйу, зует оЁ Егеапо/т П\еодга! едиайоп$ о {Те 5есопа Ктпа, ге/айуе уо/ите Гипсбоп. 





' Тпе гесеагсп {5 допе ми {Не 'ате ог {пе паерепаепе В&О. 
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